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1.  INTRODUCTION 

Prostate  cancer  (PCa)  is  the  second  leading  cause  of  cancer  death  in  men  after  lung  cancer. 
Approximately  one  in  nine  men  will  be  diagnosed  with  prostate  cancer  but  when  detected  early  and 
treated  promptly  the  five-year  relative  survival  rate  approaches  100%.  The  motivation  of  this 
project  is  to  develop  improved  means  for  detecting  prostate  cancer.  Magnetic  resonance  imaging 
(MRI)  has  been  applied  to  the  imaging  of  prostate  cancer  for  several  decades.  The  typical  MRI 
prostate  exam  today  consists  of  several  “pulse  sequences:”  (i)  T2-weighted  spin-echo  imaging;  (ii) 
diffusion- weighted  imaging  (DWI);  (iii)  dynamic  contrast-enhanced  (DCE)  perfusion  MRI.  While 
prostate  cancer  can  be  visualized  using  each  sequence,  only  sequence  (iii)  provides  dynamic 
information  about  the  temporal  enhancement  pattern  of  any  PCa  lesions.  The  purpose  of  this 
project  is  to  develop  lOx  improved  spatiotemporal  resolution  DCE-MRI  of  prostate. 


2.  KEYWORDS 

CAPR 

CE-MRA 

DCE-MRI 

MRI 

PCa 

SENSE 

SNR 


Cartesian  Acquisition  with  Projection  Reconstruction-like  sampling 
Contrast-Enhanced  Magnetic  Resonance  Angiography 
Dynamic-Constrast-Enhanced  Magnetic  Resonance  Imaging 
Magnetic  Resonance  Imaging 
Prostate  Cancer 

Sensitivity  Encoding  (a  type  of  MRI  acceleration  technique) 
Signal-to-Noise  Ratio 


3.  ACCOMPLISHMENTS: 


What  were  the  major  goals  of  the  project? 

Tasks  for  Months  1-21  (encompassing  August  15,  2015  through  June  1,  2017)  taken  from  the 
grant  application  Statement  of  Work  are  shown  below.  Completed  tasks  reported  in  last  year’s 
report  (July  2016)  are  shown  with  the  month  of  completion  in  blue  (in  parentheses).  All  tasks 
newly  completed  since  then  and  reported  here  are  shown  with  the  completion  month  shown  in 

red  (in  parentheses). 


Specific  Aim  1:  Development  of  MRI  Acquisition  Method 

Months 

Investigator 

Major  Task  1:  Optimization  of  pCAPR  Pulse  Sequence  Applied 
to  Prostate  Imaging 

Subtask  1.1:  Determine  parameter  options  for  various 
spatiotemporal  resolution  combinations 

1-6 

(5) 

Dr.  Riederer;  Dr. 
Kawashima;  Mr. 
Borisch 

Subtask  1.2:  Design,  construct,  and  test  phantom  which  mimics 
geometry  for  prostate  MRI 

1-6 

(6) 

Dr.  Riederer; 
Mr.  Hulshizer 

Subtask  1.3:  Experimentally  test  and  evaluate  versions  of  pulse 
sequence  with  prostate  phantom  and  select  optimum 
sequence  with  standard  receiver  coil  array 

6-12 

(9, ongoing) 
(12) 

Dr.  Riederer; 
Mr.  Hulshizer 

Subtask  1.4:  Apply  initial  optimized  pulse  sequence  to  three 
volunteers  with  standard  receiver  coil  array 

13-15 

(15) 

Dr.  Riederer; 

Dr.  Kawashima 

5 


Subtask  1.5:  Experimentally  test  and  evaluate  pulse  sequence  and 
triangular  array  with  prostate  phantom  and  select 
optimum  sequence 

16-24 

(20) 

Dr.  Riederer; 

Mr.  Hulshizer 

Subtask  1.6:  Determine  parameter  options  for  pulse  sequence 
using  prototype  triangular  coil  array 

18-24 

(20,  ongoing) 

Dr.  Riederer; 

Mr.  Borisch 

Specific  Aim  2:  Development  of  Special  Purpose  Receiver  Coil 
Arrays 

Major  Task  2:  Optimization  of  pCAPR  Pulse  Sequence  Applied 
to  Prostate  Imaging 

Subtask  2.1:  Select  optimum  triangular  element  size  and 
construct  matched  single  pair  of  coils 

1-6 

(3) 

Dr.  Riederer; 
Mr.  Hulshizer 

Subtask  2.2:  Construct,  tune,  and  match  multiple  pairs  of 
triangular  coils 

7-12 

(9) 

Dr.  Riederer; 
Mr.  Hulshizer 

Subtask  2.3:  Test  modular  triangular  coil  element  prototype  array  with 
up  to  16  elements  (8  coil  pairs)  in  experimental  phantom  studies 

13-18 

(18) 

Dr.  Riederer; 

Mr.  Hulshizer 

Specific  Aim  3:  Formation  of  Optimized  pCAPR  Images 

Major  Task  3.1:  Image  Reconstruction  and  System  Integration 

Subtask  3.1:  Allowance  for  arbitrary  acceleration  factors 
(Ry,Rz)  and  CAIPIRINHA  kernels  for  arbitrary  R 

1-6 

(6) 

Mr.  Borisch 

Subtask  3.2:  Incorporate  pCAPR  pulse  sequence  into  image 
reconstruction  framework 

7-12 

(7) 

Mr.  Borisch 

Subtask  3.3:  Implement  multi-processor  execution  of  individual 
modules  of  reconstruction  process 

12-24 

(18) 

Mr.  Borisch 

What  was  accomplished  under  these  goals? 

The  grant  activities  were  subdivided  into  three  specific  aims  with  subtasks  for  each.  Progress 
has  been  made  for  each  aim.  The  following  paragraphs  are  associated  with  the  specific  subtasks 
identified  in  the  Statement  of  Work  (SOW). 

Specific  Aim  1 

Subtask  1.3.  Work  on  this  subtask  continued  since  the  time  of  last  year’s  progress  report.  The 
phantom  previously  developed  was  further  modified  to  include  several  simulated  lesions  within  the 
inclusion  simulating  the  prostate.  A  sagittal  MR  image  of  the  phantom  is  shown  in  Figure  1.  The 
overall  size  of  the  phantom  simulates  a  patient  of  an  approximate  BMI  of  23.  A  larger  phantom  is 
also  being  constructed.  Related  to  this  we  have  studied  the  effect  of  temporal  resolution  on  the 
fidelity  with  which  perfusion  parameters  are  estimated.  This  was  done  by  starting  with  the  data 
acquired  at  6.5  sec  frame  time  and  down  sampling;  i.e.  taking  every  other  point  or  possibly  every 
third  point  Results  similar  to  those  in  Figure  2  show  that  temporal  sampling  coarser  than  the  6.5 
sec  used  in  this  work  causes  underestimation  of  perfusion  parameters. 

Subtask  1.4.  This  subtask  involved  use  of  the  optimized  pulse  sequence  in  three  volunteers 
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(min'1),  Case  :  1,  Slice  :  19 


Figure  2.  Plot  illustrating  (A)  downsampled  even 
and  odd  time  points  and  maps  of  Ktrans  using  all 
data  (6.5  sec)  (B)  and  only  the  odd  timepoints  (C). 

.  0  .  ,  ,  „  .  .  The  latter  markedly  causes  underestimation  of 

Figure  1.  Sagittal  MR  image  of  Ktrans 

prostate  phantom  with  simulated 

lesion.  using  the  standard  receiver  coil  array.  The  principal 

target  group  of  patients  is  that  which  is  “treatment- 
naive”  in  that  prostate  cancer  is  suspected  but  the  presence  or  degree  of  disease  must  be  determined 
prior  to  determination  of  possible  treatment.  The  optimum  sequence  was  identified  in  the  first 
year,  and  for  such  patients  temporal  sampling  the  realm  of  5  to  7  seconds  represented  a  good 
tradeoff  between  spatial  and  temporal  resolution.  Volunteers  have  been  imaged  using  the  standard 
12-element  receiver 
coil  array,  and  in  fact 
this  has  been 
expanded  to  include 
studies  of  patients 
who  are  referred  for 
prostate  MRI. 

Further,  we  have 
studied  to  possible 
extension  of  the 
acquisition  to  include 
more,  specifically  32, 
of  the  multicoil 
elements  available. 

These  additional 
elements  are 

generally  positioned 
more  lateral  than  the 
reference  12 

elements.  For  the 
same  scan,  use  of  32 
elements  consistently 
provides  improved 
image  quality. 

Sample  results  are 
shown  in  Figure  3, 
and  we  have 
published  these 


Figure  3.  Comparison  of  images  of  the  prostate  acquired  using  12- 
element  standard  coil  (left)  and  32-element  coil  (right).  Prostate  is 
identified  within  the  dashed  ellipse.  The  right  image  has  superior 
signal-to-noise  ratio. 


findings. 


Subtask  1.5.  For  this  task  we  used  a  12-element  array  composed  of  the  new  triangular  elements 


and  imaged  a  prostate  phantom  to  assess  performance  using  the  target  prostate  perfusion  pulse 
sequence.  Sample  results  are  shown  in  Figure  4.  This  figure  shows  maximum  intensity  projection 
(MIP)  images  of  the  g-factor, 
a  measure  of  the  noise 
amplification  of  the  coil  array 
for  sagittal  (A),  coronal  (B) 
and  axial  (C)  orientations. 

Also  shown  in  (D)  is  the 
cumulative  g-factor  as 
measured  over  the  3D  volume 
being  imaged.  This  work  is 
ongoing,  and  further 
measurements  and 

comparisons  will  be  made 
with  the  standard  receiver  coil. 

Specific  Aim  2 


Subtask  2.3.  The  results  in 
Figure  3  were  generated  using 
a  12-element  coil.  As 
specified  in  this  subtask,  we 
have  further  constructed  a  16- 
element  coil  composed  of  the 

triangular  elements.  The  design  using  two-element  modules  allows  the  total  number  of  elements 


Figure  4.  Plot  of  g-factors  for  12-element  coil  composed  of 
triangular  elements. 


to  be  matched  to  the  body  habitus,  i.e.  patient  diameter, 
shown  in  Figure  5. 


A  photograph  of  the  16-element  coil  is 


Specific  Aim  3 


Subtask  3.3.  The  reconstruction  hardware  was  identified  last  year  and  is  shown  in  Figure 
6.  We  have  made  progress  on  this  aim  and  have  updated  the  hardware  to  include  two  CPU 
packages,  each  with  14  physical  cores  and  28  threads,  and  thus  a  total  of  28  cores  (processor 
elements)  and  56  threads.  This  routinely  allows  fast,  clinically  acceptable  reconstruction 
times  of  3D  time-resolved,  multi-coil  (32  element)  accelerated  images  of  dynamic-contrast- 
enhanced  (DCE)  MRI.  Each  of  the  steps  in  the  reconstruction  (determination  of  coil  sensitivities, 
Fourier  transformation  of  raw  data,  un-aliasing  of  reconstructed  images,  presentation  into  3D  time 
series,  conversion  to  vendor-compatible  format)  comprises  an  individual  module. 


8 


MRI  System  Reconstruction  Hardware 

.  *?"  ISpSflH 

, - .  Parallel 


Recon 


Figure  6.  Schematic  showing  host  MRI  scanner  (left)  and  custom 
computation  hardware  (right). 


What  opportunities  for  training  and  professional  development  has  the  project  provided? 

Nothing  to  report. 

How  were  the  results  disseminated  to  communities  of  interest? 

Nothing  to  report  other  than  that  in  Section  6  of  this  report. 

What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the 
goals? 

We  will  continue  our  work  according  to  the  tasks  shown  in  the  original  Statement  of  Work. 


4.  IMPACT:  Describe  distinctive  contributions,  major  accomplishments,  innovations,  successes,  or 
any  change  in  practice  or  behavior  that  has  come  about  as  a  result  of  the  project  relative  to: 

What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

The  tasks  completed  in  the  first  21  months  of  this  three- year  project  provide  a  good  basis  for 
ongoing  work  for  all  three  Aims  in  the  months  ahead. 

What  was  the  impact  on  other  disciplines? 

Nothing  to  Report 

What  was  the  impact  on  technology  transfer? 

Our  baseline  DCE-MRI  pulse  sequences  are  now  used  clinically  at  Mayo  Clinic,  both  in  Rochester 
and  at  Mayo  Clinic  Arizona.  This  is  approximately  ten  patient  studies  per  day. 

What  was  the  impact  on  society  beyond  science  and  technology? 

Nothing  to  Report 
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5.  CHANGES/PROBLEMS:  The  PD/PI  is  reminded  that  the  recipient  organization  is  required  to 
obtain  prior  written  approval  from  the  awarding  agency  grants  official  whenever  there  are 
significant  changes  in  the  project  or  its  direction.  If  not  previously  reported  in  writing,  provide  the 
following  additional  information  or  state,  “Nothing  to  Report,”  if  applicable: 

Changes  in  approach  and  reasons  for  change 

Nothing  to  Report 

Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them 

Nothing  to  Report 

Changes  that  had  a  significant  impact  on  expenditures 

Nothing  to  Report 

Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animals,  biohazards,  and/or 
select  agents 

Significant  changes  in  use  or  care  of  human  subjects 

Nothing  to  Report 

Significant  changes  in  use  or  care  of  vertebrate  animals 

Nothing  to  Report 

Significant  changes  in  use  of  biohazards  and/or  select  agents 

Nothing  to  Report 


6.  PRODUCTS:  List  any  products  resulting  from  the  project  during  the  reporting  period.  If 
there  is  nothing  to  report  under  a  particular  item,  state  “Nothing  to  Report.” 

•  Publications,  conference  papers,  and  presentations 

Report  only  the  major  publication(s)  resulting  from  the  work  under  this  award. 

Journal  publications. 

1.  Riederer  SJ,  Borisch  EA,  Froemming  AT,  Grimm  RC,  Kawashima  A,  Mynderse  LA, 
Trzasko  JD,  “Improved  performance  of  prostate  DCE-MRI  using  a  32-coil  vs.  12-coil 
receiver  array.”  Magn  Reson  Imaging  39:15-23  (2017). 

2.  Kargar  S,  Borisch  EA,  Froemming  AT,  Kawashima  A,  Mynderse  LA,  Stinson  EG, 
Trzasko  JD,  Riederer  SJ,  “Robust  and  efficient  pharmacokinetic  parameter  estimation 
for  dynamic  contrast  enhanced  MRI  of  the  prostate,”  (in  preparation). 


Books  or  other  non-periodical,  one-time  publications. 

Nothing  to  Report 
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Other  publications,  conference  papers  and  presentations. 

1.  Kargar  S,  Stinson  EG,  Borisch  EA,  Froemming  AT,  Kawashima  A,  Mynderse  LA,  Trzasko  JD, 
Riederer  SJ,  “Patient- specific  vs.  population-based  arterial  input  function  in  perfusion  estimation 
for  DCE-MRI  of  the  prostate.”  Presented  at  28th  Annual  Meeting  of  Society  of  Magnetic 
Resonance  Angiography,  Chicago  IL,  September  2016. 

2.  Kargar  S,  Borisch  EA,  Froemming  AT,  Kawashima  A,  Mynderse  LA,  Stinson  EG, 
Trzasko  JD,  Riederer  SJ,  “Robust  and  efficient  pharmacokinetic  parameter  estimation 
for  dynamic  contrast  enhanced  MRI  of  the  prostate."  Presented  at  24th  Annual 
Meeting,  ISMRM,  Honolulu  HI,  April  2017. 

3.  Website(s)  or  other  Internet  site(s) 

Nothing  to  Report 

4.  Technologies  or  techniques 

Nothing  to  Report 

5.  Inventions,  patent  applications,  and/or  licenses 

Nothing  to  Report 

6.  Other  Products 

Nothing  to  Report 


7.  PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 


What  individuals  have  worked  on  the  project? 


Name: 

Project  Role: 

Nearest  person  month  worked: 
Contribution  to  Project: 


Name: 

Project  Role: 

Nearest  person  month  worked: 
Contribution  to  Project: 


Name: 

Project  Role: 

Nearest  person  month  worked: 
Contribution  to  Project: 


Stephen  J.  Riederer,  Ph.D. 

Principal  Investigator 

1.8 

Dr.  Riederer  directs  all  technical  aspects  of  the 
projects. 

Akira  Kawashima,  M.D.,  Ph.D. 

Co-Investigator 

.22 

Dr.  Kawashima  oversees  feasibility  testing  performed 
in  volunteers  and  provides  feedback  on  intermediate 
results  for  all  projects. 

Eric  A.  Borisch 

Information  Services  Technical  Specialist 
2.76 

Mr.  Borisch  is  responsible  for  writing  and 
developing  production-level  software  for  all 
projects,  software  heavily  centered  on 
reconstruction  of  2D-accelerated  3D  data  sets 
acquired  with  various  view  orders. 
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Name:  Thomas  C.  Hulshizer 

Project  Role:  MR  Technician 

Nearest  person  month  worked:  1 .2 

Contribution  to  Project:  Mr.  Hulshizer  is  responsible  for  construction  of  the 

pelvis-prostate  phantom,  construction  and  tuning  of 
RF  coils,  and  testing  of  prototype  MR  pulse 
sequences  using  phantoms. 


Has  there  been  a  change  in  the  active  other  support  of  the  PD/PI(s)  or  senior/key  personnel 
since  the  last  reporting  period? 

Stephen  J.  Riederer,  Ph.D. 

Ended  February  2017 
NIH  R01  EB0002 12-27 

“Modified  Elliptical  Centric  View  Orders  for  Improved  Real-Time  MRA”  (PI:  SJ  Riederer,  Ph.D.) 
Purpose:  develop  image  processing  methods  for  estimating  prostate  perfusion  parameters 
3.6  months/year  (30%  effort) 


What  other  organizations  were  involved  as  partners? 

Nothing  to  Report  other  than  the  general  support  provided  by  the  Pi’s  institution,  Mayo  Clinic. 


8.  SPECIAL  REPORTING  REQUIREMENTS 

None 


9.  APPENDICES:  Attach  all  appendices  that  contain  information  that  supplements,  clarifies  or 
supports  the  text.  Examples  include  original  copies  of  journal  articles,  reprints  of  manuscripts  and 
abstracts,  a  curriculum  vitae,  patent  applications,  study  questionnaires,  and  surveys,  etc. 

Manuscripts 

Attached  is  a  pdf  copy  of  the  published  article: 

“Improved  Performance  of  Prostate  DCE-MRI  Using  a  32-Coil  vs.  12-Coil  Receiver  Array” 


Magnetic  Resonance  Imaging  39  (2017)  15-23 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Magnetic  Resonance  Imaging 


journal  homepage:  www.mrijournal.com 


Original  contribution 
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Purpose:  To  assess  whether  acquisition  with  32  receiver  coils  rather  than  the  vendor-recommended  12  coils  pro¬ 
vides  significantly  improved  performance  in  3D  dynamic  contrast-enhanced  MRI  (DCE-MRI)  of  the  prostate. 
Materials:  The  study  was  approved  by  the  institutional  review  board  and  was  compliant  with  HIPAA.  50  consec¬ 
utive  male  patients  in  whom  prostate  MRI  was  clinically  indicated  were  prospectively  imaged  in  March  2015 
with  an  accelerated  DCE-MRI  sequence  in  which  image  reconstruction  was  performed  using  12  and  32  coil 
elements.  The  two  reconstructions  were  compared  quantitatively  and  qualitatively.  The  first  was  done  using  sig- 
nal-to-noise  ratio  (SNR)  and  g-factor  analysis  to  assess  sensitivity  to  acceleration.  The  second  was  done  using  a 
five-point  scale  by  two  experienced  radiologists  using  criteria  of  perceived  SNR,  artifact,  sharpness,  and  overall 
preference.  Significance  was  assessed  with  the  Wilcoxon  signed  rank  test.  Extension  to  T2-weighted  spin-echo 
and  diffusion  sequences  was  assessed  in  phantom  studies. 

Results:  Reconstruction  using  32  vs.  12  coil  elements  provided  improved  performance  in  DCE-MRI  based  on 
intrinsic  SNR  (18%  higher)  and  g-factor  statistics  (14%  higher),  with  a  median  32%  higher  overall  SNR  within 
the  prostate  volume  over  all  subjects.  Reconstruction  using  32  coils  was  qualitatively  rated  significantly  im¬ 
proved  (p  <  0.001 )  vs.  1 2  coils  on  the  basis  of  perceived  SNR  and  radiologist  preference  and  equivalent  for  sharp¬ 
ness  and  artifact.  Phantom  studies  suggested  the  improvement  in  intrinsic  SNR  could  extend  to  T2-weighted 
spin-echo  and  diffusion  sequences. 

Conclusions:  Reconstruction  of  3D  accelerated  DCE-MRI  studies  of  the  prostate  using  32  independent  receiver 
coils  provides  improved  overall  performance  vs.  using  12  coils. 

©  2017  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

MR  imaging  of  the  prostate  is  commonly  performed  using  a  multi- 
parametric  approach  in  which  multiple  sequences  are  used  to  aid  in  ra¬ 
diologic  interpretation  [1-3].  In  addition  to  T2-weighted  spin-echo 
(T2SE)  and  diffusion-weighted  image  (DWI)  a  pulse  sequence  typically 
used  within  this  exam  is  three-dimensional  (3D)  dynamic  contrast- 
enhanced  MRI  (DCE-MRI)  [4-6]  in  which  a  contrast  agent  is  adminis¬ 
tered  intravenously,  and  images  are  acquired  of  the  prostate  to  observe 
washin  and  washout  of  the  contrast-enhanced  blood  over  the  entire 
prostate  volume.  DCE-MRI  continues  to  be  used  in  the  context  of  the 
PI-RADS  version  2  reporting  system  for  prostate  cancer  for  treatment- 
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naive  prostate  glands  [7[.  Further,  recurrence  of  prostate  cancer  after 
definitive  treatment  has  been  reported  to  be  in  the  15  to  40%  range  [8, 
9],  and  in  such  cases  DCE-MRI  has  been  shown  to  be  of  high  value  [10, 
11].  Because  the  desired  spatiotemporal  resolution  of  the  DCE-MRI  se¬ 
quence  typically  pushes  the  limits  of  signal-to-noise  ratio  (SNR),  the 
need  for  good  performance  of  the  receiver  coils  is  important  [12  j. 

For  prostate  MRI  the  coil  usage  reported  in  the  literature  is  highly 
variable.  One  significant  choice  is  whether  an  endorectal  (ER)  coil  should 
be  used,  and  several  studies  have  compared  ER  with  non-ER  acquisition 
[13-16].  However,  whether  an  ER  coil  is  used  or  not,  in  contemporary 
prostate  MRI  some  kind  of  multi-element  receiver  coil  is  generally  used. 
A  sampling  of  the  recent  literature  shows  use  of  multi-element  coils  hav¬ 
ing  four  [17],  six  [18],  eight  [11,15,19],  12  to  15  [6],  and  18  [20]  elements, 
or  with  the  number  of  elements  not  explicitly  given  [4,5], 

For  many  contemporary  MRI  systems  the  number  of  receiver  channels 
available  is  32  or  more.  Also,  coil  arrays  with  literally  dozens  of  elements 
are  now  embedded  within  the  patient  table.  Although  both  of  these  fac¬ 
tors  can  greatly  facilitate  the  usage  of  a  large  number  of  elements  for  a 


http://dx.doi.org/10.1016/jjnri2017.01.017 
0730-725X/©  2017  Elsevier  Inc.  All  rights  reserved. 


16 


S.J.  Riederer  et  al.  /  Magnetic  Resonance  Imaging  39  (2017)  15-23 


given  exam,  it  is  not  clear  what  the  number  or  combination  of  elements 
should  be  for  prostate  MRI.  If  an  ER  coil  is  not  used,  then  it  is  important 
to  have  a  multicoil  setup  which  provides  high  SNR  in  the  region  of  the 
prostate.  The  results  of  Ref.  [16]  indicated  superior  performance  for  the 
combination  of  an  ER  coil  and  a  16-element  multicoil  vs.  use  of  a  six-ele¬ 
ment  multicoil  alone.  However,  separation  of  the  improvement  due  to  the 
change  of  the  multicoil  vs.  use  of  the  ER  coil  was  not  determined. 

Prostate  MRI  at  our  institution  is  principally  performed  using  a  wide- 
bore  3  T  MRI  system  with  features  similar  to  those  described  above.  The 
vendor-recommended,  default  coil  selection  for  prostate  DCE-MR1  calls 
for  12  channels  of  data  acquisition  to  be  used  of  the  32  available  on  the 
system.  These  12  channels  include  eight  active  elements  from  the  array 
contained  within  the  table  posterior  to  the  supine  patient  and  four  ac¬ 
tive  elements  used  from  a  16-element  array  placed  anteriorly.  In  inves¬ 
tigating  to  what  extent  the  DCE-MRI  sequence  could  be  accelerated,  we 
considered  whether  additional  coil  elements  could  be  used.  To  fully  ex¬ 
ploit  the  capability  of  the  MRI  system,  the  option  for  use  of  all  32  chan¬ 
nels  was  considered.  Thus,  the  specific  hypothesis  of  this  work  was  that 
the  use  of  32  vs.  1 2  independent  receiver  coils  would  provide  improved 
performance  in  DCE-MRI  of  the  prostate.  The  feasibility  of  improved 
performance  for  other  sequences  of  a  multi-parametric  prostate  MRI 
exam  was  also  assessed  with  phantom  studies. 

2.  Materials  and  methods 

This  study  was  approved  by  the  institutional  review  board  which 
waived  the  need  for  written  consent.  The  study  was  compliant  with 
HIPAA. 

2.1.  Subjects 

50  consecutive  male  subjects  for  whom  a  prostate  MRI  exam  was 
clinically  indicated  and  who  gave  their  assent  for  their  exam  results  to 
be  used  for  research  purposes  were  prospectively  enrolled  in  the 
study  over  the  period  March  10-27,  2015.  The  age,  weight,  and  body 
mass  index  (BMI)  ranges  were  51  to  86  years,  63.5  to  155.1  kg,  and 
22.5  to  46.8,  respectively.  Forty-four  of  the  50  had  intact  prostates;  six 
were  evaluated  post-prostatectomy.  Other  information  on  the  patient 
cohort,  including  radiological  interpretation,  is  shown  in  Table  1.  Of 
note  is  that  34%  ( 1 7/50)  of  the  patients  were  imaged  as  followup  to  pre¬ 
vious  treatment  for  prostate  cancer. 

2.2.  MRI  acquisition 

All  studies  were  performed  on  either  of  two  identical  3.0  T  MRI  scan¬ 
ners  (Discovery  MR750w,  GE  Healthcare,  Waukesha  WI)  utilizing  an  in¬ 
stitutional  clinical  exam  protocol.  Each  machine  has  a  70  cm  diameter 
bore,  a  vendor-provided  40-element  receiver  coil  array  (Geometry 
Embracing  Method  “GEM"  array)  embedded  within  the  patient  table, 
and  32  receiver  channels.  Each  patient  exam  included  a  localizer,  T2- 
weighted  spin-echo,  and  diffusion-weighted  sequences,  followed  by  a 
DCE-MRI  study  performed  with  intravenous  contrast  administration. 


Details  of  the  RF-spoiled  gradient  echo  DCE-MRI  sequence  are  shown 
in  Table  2.  Contrast  material  (Dotarem,  Guerbet,  Paris,  France)  was  ad¬ 
ministered  into  an  arm  vein  at  a  rate  of  3  ml/s  followed  by  a  20  ml  saline 
flush  at  3  ml/s.  The  contrast  dose  was  0.1  mmol/kg,  with  a  maximum  of 
20  ml  for  patients  weighing  100  kg  or  more. 

The  impetus  for  this  study  was  to  push  the  acceleration  of  the  DCE- 
MRI  sequence  for  improved  spatiotemporal  resolution.  The  sequence 
used  is  based  on  one  developed  for  time-resolved  contrast-enhanced 
MR  angiography  (CE-MRA)  using  two-dimensional  (2D)  SENSE  acceler¬ 
ation  [21  ]  and  view  sharing  [22].  For  this  work  the  slab  orientation  was 
approximately  axial  but  with  slight  forward  tilting  to  align  the  slab  se¬ 
lect  direction  with  the  central  axis  of  the  prostate  gland  as  determined 
in  the  sagittal  localizer.  SENSE  acceleration  factors  of  2.49  and  1.12 
were  applied  along  the  left/right  (L/R)  phase  encode  and  approximate 
superior/inferior  (S/I)  slice  encode  directions,  respectively,  yielding  a 
net  acceleration  factor  of  R  =  2.78.  DCE-MRI  acquisition  was  initiated 
20  s  prior  to  the  start  of  contrast  injection,  the  frame  time  was  approx¬ 
imately  6.6  s,  and  a  total  of  33  time  frames  were  collected.  Versions  of 
the  sequence  with  higher  accelerations  were  considered  for  shorter 
frame  times;  e.g.  R  =  5.0  for  4.0  s  frame  time,  but  for  the  same  spatial 
resolution  the  SNR  loss  was  considered  too  severe. 

2.3.  Selection  of  receiver  coil  elements 

This  work  made  use  of  the  40-element  GEM  receiver  coil  array 
shown  schematically  in  Fig.  1A  with  a  male  torso  and  prostate  shown 
approximately  to  scale.  The  array  consists  of  five  columns  of  elements 
oriented  longitudinally  within  the  patient  table  which  for  the  supine  pa¬ 
tient  are  located  posteriorly.  Also  used  was  a  16-element  coil  array 
placed  anteriorly  consisting  of  four  longitudinally  oriented  columns 
each  comprised  of  four  elements  (Fig.  IB).  In  both  figures  an  approxi¬ 
mate  75  cm3  sphere,  simulating  the  prostate  gland,  is  shown  in  red  to 
give  a  sense  of  scale  of  typical  gland  size  to  extent  of  coil  coverage. 

The  combinations  of  coil  elements  available  for  usage  are  limited  by 
the  vendor,  and  recommendations  are  made  according  to  the  type  of 
exam.  For  prostate  MRI  this  calls  for  12  active  receiver  coil  elements, 
eight  from  the  GEM  array  and  four  from  the  anterior  array.  These  are 
highlighted  in  yellow  in  Figs.  1A  and  B,  respectively.  The  elements  se¬ 
lected  from  the  GEM  array  (A)  are  the  two  central-most  elements 
which  encompass  the  S/I  extent  of  the  prostate,  as  identified  from  the 
sagittal  scout  images,  and  the  three  elements  from  the  next  closest  col¬ 
umns  on  each  side.  For  the  anterior  coil  (B)  the  two  central  elements 
from  the  two  rows  which  similarly  encompass  the  S/I  extent  of  the  pros¬ 
tate  are  selected.  All  other  coil  elements  from  both  arrays  are  electroni¬ 
cally  disabled  during  acquisition,  and  data  from  the  12  elements  are 
individually  digitized  and  used  in  reconstruction. 

To  attempt  to  exploit  the  full  32-channel  capability  of  the  MRI 
systems  for  accelerated  DCE-MRI,  we  next  considered  use  of  the  ven¬ 
dor-allowed  1 6  elements  from  each  array.  This  is  also  depicted  in  Figs. 
1A-B.  Coil  elements  used  for  the  12-element  case  were  expanded  to  in¬ 
clude  those  shown  in  blue.  For  the  GEM  coil  (A)  the  two  central-most  el¬ 
ements  best  aligned  with  the  S/I  prostate  extent  are  selected  as  before. 


Table  1 

Summary  of  patient  information  and  radiological  interpretation.  Interpretation  was  based  on  the  full  multi-parametric  MRI  exam,  not  just  the  DCE-MRI  sequence. 


Patients  with  no  prior  treatment 

Number 

Interpretation: 

Interpretation: 

no  suspicion  of  significant  prostate  cancer 

suspicion  of  significant  prostate  cancer 

33 

16 

17 

Patients  with  prior  treatment 

Number 

Interpretation: 

Interpretation: 

no  worrisome 

worrisome 

lesion 

lesion 

Radical  prostatectomy 

6 

2 

4 

Radiation  therapy  (external  beam  or  brachytherapy, 

10 

2 

8 

including  with  possible  adjuvant  hormonal  therapy) 

Hormonal  therapy  only 

1 

0 

1 

Total  (patients  with  prior  treatment) 

17 

4 

13 

Total  (all  patients) 

50 

20 

30 
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Table  2 

Parameters  for  3D  RF-spoiled  gradient  echo  prostate  DCE-MRI  sequence.  Acquisition 
of  coil  calibration  images  used  the  same  sequence  applied  once  with  phase  resolution 
reduced  from  384  to  192. 


Parameter 

Value 

Repetition  time  (TR) 

5.3  msec 

Echo  time  (TE) 

2.2  msec 

aField  of  view 

220  x  440  x  114  mm3 

Sampling  resolution 

256  x  384  x  38 

aSpatial  resolution 

0.86  x  1.15  x  3.0  mm3 

Acceleration 

2.49  (Ry)  x  1.12 

(Rz)  =  2.78 

Frame  time 

6.6  s 

Temporal  footprint 

^19  s 

Number  of  frames 

33 

Scan  time 

Rs  3.5  min 

a  These  parameters  are  all  expressed  as  (frequency  x  phase  x  slice)  =  (A/P  xI/Rx  S/I). 


The  lateral  rows  containing  those  elements  are  supplemented  with  the 
rows  of  elements  next  positioned  superiorly  and  inferiorly.  For  each  row 
the  leftmost  and  rightmost  elements  (e.g.  elements  6  and  36  of  Fig.  1A) 
are  automatically  combined  in  vendor  hardware  prior  to  digitization,  in 
effect  forming  one  virtual  coil  from  two  coil  elements.  Thus,  the  20  ele¬ 
ments  contained  within  the  four  selected  rows  of  the  GEM  coil  are 
encoded  in  16  individual  coils,  with  four  of  these  being  two-element 
combinations.  For  the  anterior  array  all  16  elements  are  used  (B).  It  is 
noted  that  the  channels  used  for  the  12-coil  acquisition  are  a  subset  of 
those  used  for  the  32-coil  reconstruction.  The  remaining  20  coil  ele¬ 
ments  located  at  the  ends  of  the  GEM  array  are  electronically  disabled 
during  the  32-coil  acquisition. 

2.4.  Image  reconstruction 

Reconstruction  was  performed  offline  with  standard  SENSE 
unfolding  [23]  using  a  custom-built  computing  system  described  in 
[24].  To  avoid  the  complications  and  potential  variability  of  performing 
two  separate  DCE-MRI  studies  on  each  subject,  we  investigated  if  a  sin¬ 
gle  32-coil  acquisition  could  be  done  and  the  12-coil  acquisition  accu¬ 
rately  simulated  by  using  only  the  appropriate  12  data  sets  for 
reconstruction.  The  risk  with  this  approach  is  that  the  electronically  ac¬ 
tive  but  unused  20  coil  elements  in  the  32-active-element  acquisition 
would  interfere  through  undesirable  coupling  with  the  12  elements  se¬ 
lected  for  reconstruction.  To  assess  this  we  first  performed  test  scans  in 
a  volunteer  in  which  separate  acquisitions  were  done  with  the  12-ac- 
tive-coil  and  32-active-coil  approaches  using  the  accelerated  DCE-MRI 
sequence  without  contrast  injection.  Unaccelerated  coil  calibration 
image  sets  were  also  acquired  with  both  approaches.  Data  from  the 
32-coil  acquisition  were  reconstructed  two  ways:  (i)  using  all  32  coils; 
and  (ii)  using  data  from  only  the  same  12  coils  as  for  the  12-active- 
coil  acquisition.  Reconstruction  (iii)  was  done  using  all  12  coils  of  the 
12-active-coil  acquisition. 

The  images  from  reconstructions  (i),  (ii),  and  (iii)  were  compared  in 
two  ways.  First,  images  of  absolute  SNR  were  formed  from  data  sets  (i), 
(ii),  and  (iii)  using  the  method  of  Refs  [25,26].  As  the  coil  sensitivity  pro¬ 
files  used  in  these  calculations  were  estimated  empirically  via  the  root- 
sum-of-square  demodulation,  these  SNR  values  are  quantitatively  ap¬ 
proximate.  Reconstructed  SNR  values  were  taken  of  a  3D  volume  just 
encompassing  the  prostate,  in  this  case  approximately  65  cm3  and  com¬ 
prised  of  >20,000  pixels.  Second,  using  the  coil  calibration  data  and  as¬ 
suming  the  acceleration  factors  of  the  DCE-MRI  sequence  in  Table  2, 
images  were  formed  of  the  g-factor,  a  mathematical  measure  of  the  abil¬ 
ity  of  a  receiver  coil  array  to  retain  SNR  in  accelerated  MR  acquisition 
[23],  These  comparisons,  shown  in  Fig.  2,  indicated  that  the  SNR  (A) 
and  g-factor  statistics  (B)  of  reconstructions  (ii)  and  (iii)  were  essential¬ 
ly  indistinguishable,  and  both  were  different  from  results  for  recon¬ 
struction  (i).  Consequently,  results  for  the  12-coil  reconstruction  were 
generated  in  the  patient  study  by  selecting  data  only  from  those  12 


elements  from  the  32-active-coil  acquisition  and  reconstructing  that 
data.  Comparisons  were  then  made  with  the  reconstruction  using  the 
full  32-coil  data  set. 

2.5.  Radiological  evaluation 

For  each  of  the  50  patient  studies  the  32-coil  DCE-MRI  sequence  was 
reconstructed,  and  images  were  zero  padded  by  a  factor  of  2  and 
cropped  to  effectively  magnify  the  region  surrounding  the  prostate.  An 
observer  not  performing  the  radiological  evaluation  selected  an  axial 
partition  midway  through  the  S/I  extent  of  the  prostate  at  the  time 
frame  closest  to  50  s  post  injection.  This  typically  corresponded  to  a 
time  10  to  20  s  after  peak  contrast  enhancement  of  any  rapidly  enhanc¬ 
ing  lesions  in  the  prostate.  The  magnified  32-coil  and  12-coil  images 
from  this  partition  and  time  frame  were  then  placed  side-by-side  ran¬ 
domly  and  in  blinded  fashion  for  each  study.  This  set  of  50  composite 
images  was  then  provided  to  each  reviewer. 

The  two  radiologist  reviewers  (ATF,  six  years'  experience  in  prostate 
MRI;  AK,  20  years'  experience)  then  independently  graded  each  image 
pair  using  a  five-point  scale  (  — 2  =  left  (L)  image  significantly  better 
than  right  (R)  image;  —  1  =  L  slightly  better  than  R;  0  =  L  and  R  images 
equivalent;  + 1  =  R  slightly  better  than  L;  +  2  =  R  significantly  better 
than  L).  This  was  done  for  each  of  the  four  criteria  of  perceived  SNR, 
level  of  artifact,  sharpness,  and  overall  preference. 

In  addition,  images  of  absolute  SNR  were  made  for  the  32-coil  and 
1 2-coil  reconstructions  for  each  of  the  studies  using  the  coil  calibration 
images  as  described  previously  for  the  volunteer  study.  For  this  analysis 
images  from  the  six  of  the  50  patients  who  were  imaged  post-prostatec- 
tomy  were  excluded.  For  each  study  the  3D  rectangular  volume  was 
identified  on  the  reconstructed  images  which  just  encompassed  the 
prostate,  histograms  of  the  SNR  values  were  generated  for  all  voxels 
within  the  volume  for  the  two  reconstructions,  and  the  ratio  of  median 
values  of  the  32-coil  vs.  12-coil  histograms  was  taken  as  a  measure  of 
SNR  improvement  for  that  study.  These  volumes  ranged  from  28  to 
265  cm3  with  a  median  of  101  cm3,  typically  including  several  tens  of 
thousands  of  pixels.  This  process  was  repeated  for  the  g-factor  for  the 
acceleration  used  in  the  DCE-MRI  run,  and  the  ratio  of  median  g-factor 
values  determined. 

2.6.  Statistical  significance 

For  the  radiological  evaluation  after  accounting  for  the  blinded  pre¬ 
sentation,  significant  (defined  as  p<  0.05)  difference  from  the  null  hy¬ 
pothesis  of  equivalent  performance  was  assessed  with  the  Wilcoxon 
signed  rank  test. 

2.7.  Extension  to  other  pulse  sequences  for  prostate  MRI 

After  having  demonstrated  improved  performance  in  prostate  DCE- 
MRI,  it  was  of  interest  to  assess  potential  improvement  in  the  other 
sequences  of  multi-parametric  prostate  MRI,  namely  T2-weighted 
spin-echo  (T2SE)  and  diffusion-weighted  imaging  (DWI).  This  was  not 
possible  with  the  50  patient  studies  described  previously.  To  study 
this  in  a  controlled  manner  a  phantom  was  used,  consisting  of  a  plastic, 
male-like  pelvic  shell  filled  with  B-gel,  sorbic  acid,  and  salt  dissolved  in 
water  which  upon  congealing  mimics  soft  tissue.  During  the  pouring 
two  inclusions  were  added,  one  a  50  cm3  balloon  to  mimic  the  bladder 
and  the  second  a  40  cm3  balloon  to  mimic  the  prostate.  For  both 
balloons  the  B-gel  concentration  was  adjusted  to  provide  slight  but 
discernible  signal  differences  vs.  background.  A  photograph  of  the  phan¬ 
tom  is  shown  in  Fig.  8A.  Based  on  phantom  size,  it  corresponds  to  a  BMI 
of  approximately  24. 

The  phantom  was  imaged  with  the  12-  and  32-coil  arrangements 
using  the  sequences:  (i)  DCE-MRI  described  above,  (ii)  coronal  T2SE, 
(iii)  full  lateral-FOV  “conventional”  axial  DWI,  and  (iv)  limited-lateral- 
FOV  axial  DWI  [27],  Sequence  parameters  are  shown  in  Table  3.  The 
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(A) 


Posterior  Array  Element  Selection 


Superior  < - >  Inferior 


12-coil  operation 
32-coil  operation 


n 


n 


(B) 


Anterior  Array  Element  Selection 


»  Inferior 


□  12-coil  operation 

□  32-coil  operation 


Fig-1  .  Schematics  of  coil  elements  and  element  selection  for  1 2-coil  and  32-coil  operation.  (A)  40-element  GEM  array  contained  within  the  patient  table  and  located  posterior  to  the  supine 
patient.  (B)  Schematic  of  16-element  array  placed  anteriorly  to  the  supine  patient.  For  12-channel  operation  the  coil  elements  shown  in  yellow  are  selected,  eight  from  the  posterior  array 
and  four  from  the  anterior  array.  For  32-channel  operation  all  coils  in  yellow  and  additionally  those  in  blue  are  selected.  For  the  posterior  GEM  array  signals  from  the  lateral-most  elements 
are  paired  (4  with  34,  5  with  35,  etc.)  and  combined  in  hardware  into  one  virtual  coil  per  pair.  For  32-channel  operation  all  16  elements  of  the  anterior  array  are  used. 


DCE-MRI  data  was  used  to  reconstruct  3D  images  of  absolute  SNR  as  de¬ 
scribed  previously,  and  the  ratio  between  the  32-  vs.  12-coil  images 
taken  as  a  measure  of  absolute  intrinsic  SNR  improvement.  Each  of  the 
other  sequences  was  run  twice,  the  same  gain  settings  for  both  runs. 
For  each  sequence  an  image  from  one  of  the  runs  was  selected  of  a  sec¬ 
tion  midway  through  the  prostate  inclusion.  The  signal  S  was  taken  as 
the  mean  over  a  circular,  approximate  2  cm2  region-of-interest  (ROI) 
within  the  inclusion.  Noise  a  was  assessed  by  forming  the  difference 
image  between  the  two  acquired  images  of  the  same  anatomic  section, 


measuring  the  standard  deviation  in  the  same  ROI  as  used  for  signal,  and 
dividing  the  result  by  i/2  to  account  for  the  subtraction.  SNR  was  then 
defined  as  S/a,  and  the  ratio  of  SNRs  of  32-  and  12-coil  acquisition 
determined. 

3.  Results 

Fig.  2B  shows  the  cumulative  histograms  of  the  g-factor  values  for 
the  32-coil  and  12-coil  reconstructions.  The  solid  green,  blue,  and  red 
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32  12  (masked)  12  (native) 


(B) 


Plot  of  Cumulative  Histograms  of  g-Factor  Values 


1  1.05  1.1  1.15  1.2  1.25  1.3  1.35  1.4  1.45  1.5 

g-factor  limit 


Fig.  2.  (A)  Box  and  whisker  plots  of  the  reconstructed  SNR  values  from  the  non-contrast- 
enhanced  test  case  evaluating  reconstructions  (i)  32-coil  acquisition  with  32-coil 
reconstruction,  (ii)  32-coil  acquisition  with  12-coil  reconstruction  (“masked”),  and  (iii) 
12-coil  acquisition  with  12-coil  reconstruction  (“native”).  Each  figure  shows  the  median, 
±25%  values  (box  boundaries),  and  ±45%  boundaries  (whiskers).  Median  values  are 
28.7,  23.9,  and  23.5  (a.u.).  (B)  Plot  of  the  cumulative  g-factor  statistics  for  the  32- 
channel  acquisition  with  32-coil  (green  lines  and  curves)  and  12-coil  (blue  lines  and 
curves)  reconstructions  and  for  the  12-channel  acquisition  with  12-coil  reconstruction 
(solid  red  line).  The  three  solid  curves  are  for  the  non-contrast-enhanced  test  scans.  The 
shaded  green  and  blue  zones  show  the  ranges  of  g-factor  values  measured  across  the 
32-coil  and  12-coil  reconstructions  for  all  50  patient  studies,  respectively.  The  dotted 
green  and  blue  lines  correspond  to  the  median  values,  the  dark  shaded  zones  to  ±  12.5% 
about  the  median  and  the  light  shaded  zones  to  ±37.5%  about  the  median.  For  all  plots 
the  statistics  are  measured  over  the  3D  volume  encompassing  the  prostate. 

lines  are  for  the  unenhanced  test  scans  of  the  volunteer.  The  green  shad¬ 
ed  region  shows  the  cumulative  g-factor  histogram  for  the  32-coil  re¬ 
construction  based  upon  all  50  patient  studies.  The  dashed  green  line 
is  the  median  value;  the  dark  green  zone  encompasses  the  central 
25%,  and  the  light  green  zone  the  central  75%  at  that  cumulative  per¬ 
centage.  Similarly,  the  blue  shaded  region  shows  analogous  results  for 
the  12-coil  reconstruction.  The  more  rapid  approach  to  100%  of  the 
green  32-coil  reconstruction  indicates  the  overall  smaller  g-factor 
values  and  better  retention  of  SNR  vs.  the  two  12-coil  reconstructions 
(blue  curves  and  blue  regions  and  red  curve).  This  distinction  between 
curves  is  maintained  across  all  50  patient  studies  in  that  the  shaded 
green  and  blue  regions  are  well  separated.  The  close  match  of  the  blue 
regions  and  red  curve  indicates  the  equivalence  in  performance  of  re¬ 
construction  (ii)  and  (iii). 

Fig.  3  shows  results  from  the  qualitative  evaluation  of  the  two  radi¬ 
ologists.  32-coil  reconstruction  (positive  scores)  was  evaluated  as 


Table  3 

Technical  parameters  of  T2SE  and  DWI  pulse  sequences  used  for  phantom-based  SNR 
analysis. 


Parameter 

T2SE 

Conventional 

DWI 

Limited-FOV 

DWI 

Format 

Coronal 

Axial 

Axial 

Field  of  View  (mm2) 

220  x  220 

360  x  360 

260  x  130 

(FOVx  x  FOVy) 

Acquisition  matrix 

320  x  320 

192  x  192 

144  x  72 

(Xx  Y) 

Inplane  resolution  (mm2) 

0.688  x  0.688 

1.875  x  1.875 

1.8  x  1.8 

(Xx  Y) 

Slice  Thickness/spacing 

3/3.3 

3/0 

4/0 

(mm  for  both) 

TR  (ms) 

4951 

5888 

3925 

TE  (ms) 

96 

70 

75 

Number  of  slices 

26 

23 

16 

b-Values  (s/mm2) 

- 

100/1000 

100/1000 

Number  of  averages 

1 

1 

1 

The  T2SE  sequence  used  a  fast-spin-echo  readout  with  an  echo  train  length  of  4. 


significantly  superior  (p  <  0.001)  to  12-coil  reconstruction  using  the 
criteria  of  SNR  (A)  and  overall  preference  ( D)  by  both  reviewers  individ¬ 
ually  and  in  aggregate.  32-coil  and  12-coil  reconstructions  were  evalu¬ 
ated  as  equivalent  for  the  criteria  of  level  of  artifact  (B)  and  sharpness 
(C)  by  both  reviewers  individually  and  in  aggregate.  For  perceived 
SNR  the  two  reviewers'  scores  matched  in  38/50  cases  and  were  within 
one  value  on  the  five-point  —  2  to  +2  scale  in  all  50/50.  For  artifact 
these  corresponding  results  were  34/50  and  50/50;  for  sharpness  27/ 
50  and  46/50;  for  overall  preference  22/50  and  47/50. 

Fig.  4  shows  the  percent  improvement  in  SNR  values  provided  by  the 
32-coil  vs.  12-coil  reconstruction  plotted  vs.  the  BMI  of  the  patient.  The 
red  circles  show  the  case  with  no  acceleration,  as  determined  from  the 
images  of  SNR  using  the  coil  calibration  data.  The  median  increase 
(18%)  over  the  44  studies  is  noted  with  the  red  hashmark  on  the  ordi¬ 
nate  axis,  and  the  trend  line,  determined  by  least  squares  regression, 
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Reviewer  preference 
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20 
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Overall  preference 
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!■ 


-2-10  1  2 

Reviewer  preference 


Fig.  3.  Histograms  showing  the  results  of  the  radiological  review  for  perceived  SNR  (A), 
level  of  artifact  (B),  sharpness  (C),  and  overall  preference  (D).  For  SNR  (A)  and  overall 
preference  (D)  the  preference  for  the  32-coil  reconstruction  (positive  scores)  was 
significant  (p  <  0.001 )  for  both  reviewers  individually  and  in  aggregate.  For  artifact  (B) 
and  sharpness  (C)  there  was  no  significant  preference. 
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Plot  of  SNR  Improvement  in  32-  vs.  12-Coil  Reconstructions 


Fig.  4.  Plot  of  the  ratio  of  median  reconstructed  SNR  values  without  (red  circles)  and  with  (black  triangles)  the  additional  effect  of  g-factor  improvements  for  the  32-coil  and  12-coil 
reconstructions  plotted  vs.  BM1  of  the  patient.  In  each  case  values  were  computed  from  a  volume  encompassing  the  prostate.  The  median  increase  for  each  is  shown  in  the 
corresponding  colored  hashmark  on  the  ordinate,  and  trend  lines  of  each  with  BM1,  as  determined  from  least  square  regression,  are  also  noted. 


shown  in  red.  The  triangles  show  the  case  when  acceleration  is  addi¬ 
tionally  used  at  the  acceleration  factors  employed  in  this  study  with 
analogous  median  (32%)  and  trend  lines  shown  in  black.  The  median 
BMI  (29.1 )  for  the  50  studies  is  also  noted  as  a  hashmark  on  the  abscissa. 


Figs.  5-7  show  sample  results  from  three  patient  studies  in  which  the 
32-coil  and  12-coil  reconstructions  are  compared.  Video  VI  compares 
all  31  time  frames  for  12-  vs.  32-coil  reconstructions  for  the  study  pre¬ 
sented  in  Fig.  5. 


Fig.  5.  Comparison  of  12-coil  (left)  and  32-coil  (right)  reconstructions  of  prostate  DCE-MR1  in  a  patient  with  BMI  26.6.  Prostate  is  identified  within  the  white  ellipse.  Radiologists  #1  and  #2 
both  assigned  scores  of  (+2, 0,  0,  +2)  for  (perceived  SNR,  artifact  level,  sharpness,  overall  preference)  where  positive  scores  reflect  preference  for  the  32-coil  result. 
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Fig.  6.  Comparison  of  12-coil  (left)  and  32-coil  (right)  reconstructions  of  prostate  DCE-MRI  in  a  patient  with  BM1 30.4.  Radiologists  #1  and  #2  assigned  scores  of  (+1,0, 0,  +  1)  and  (+1,0, 
0,  +2). 


Table  4  summarizes  results  for  relative  SNR  for  the  phantom  studies 
which  show  an  SNR  increase  ranging  from  9.7%  for  T2SE  to  34.3%  for  the 
b-1000  conventional  DWI.  Fig.  8B  is  a  axial  image  of  the  phantom  show¬ 
ing  the  SNR  increase  determined  from  the  DCE-MRI  coil  calibration  data. 


The  position  of  the  prostate  inclusion  is  designated  by  the  dashed  circle. 
Fig.  8C  is  a  side-by-side  comparison  of  unsubtracted  coronal  T2SE  im¬ 
ages  of  the  phantom  using  12-  (left)  and  32-coil  (right)  acquisition  illus¬ 
trating  the  size  and  placement  of  the  bladder  and  prostate  inclusions. 


Fig.  7.  Comparison  of  12-coil  (left)  and  32-coil  (right)  reconstructions  of  prostate  DCE-MRI  in  a  patient  with  BM1 33.4  and  with  implanted  seeds  for  brachytherapy  (black  dropouts,  e.g. 
short  white  arrows).  In  this  exam  an  endorectal  coil  used  for  sequences  other  than  DCE-MRI  was  applied  within  a  gel-filled  insert  (long  white  arrow)  but  not  active  for  the  DCE-MRI 
sequence.  Radiologists  #1  and  #2  assigned  scores  of  (+1, 0, 0,  +1)  and  (0, 0,  0, 0). 
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Table  4 

Relative  SNR  of  32-  vs.  12-coil  acquisition  as  measured  in  prostate  phantom. 


Comparison 

Sequence 

Relative  SNR 

1 

T2SE 

1.097 

2 

Conventional  DW1 

(b  =  100  s/mm2, 1  average) 

1.163 

3 

Conventional  DW1 

(b  =  1000  s/mm2, 1  average) 

1.346 

4 

Limited-FOV  DWI 

(b  =  100  s/mm2, 1  average) 

1.116 

5 

Limited-FOV  DWI 

(b  =  1000  s/mm2, 1  average) 

1.338 

(A) 


Approximate  Location 
of  Bladder  S/I  Center 


4.  Discussion 

For  the  same  level  of  acceleration,  the  acquisition  and  reconstruction 
of  data  from  32  receiver  coils  encompassing  the  pelvis  provides  superior 
measured  and  perceived  SNR  in  dynamic-contrast-enhanced  prostate 
MR1  vs.  use  of  the  vendor-recommended  12  receiver  coils.  The  median 
improvement  in  overall  measured  SNR  was  32%  over  the  volume  of 
the  prostate.  Although  the  placement  of  the  1 2  receiver  elements  used 
already  encompasses  the  full  lateral  and  superior/inferior  extent  of  the 
prostate,  the  incorporation  of  data  from  additional  coils  can  provide  im¬ 
proved  performance  without  generation  of  artifact  related  to  the  in¬ 
crease  in  number  of  elements. 

There  are  possible  disadvantages  in  the  use  of  additional  receiver 
coils  in  image  reconstruction.  When  done  on  the  vendor  system,  the 
reconstruction  time  for  a  32-coil  DCE-MRI  run  is  50  s  vs.  1 5  s  on  our  cus¬ 
tom  hardware.  Lacking  high  speed  computational  hardware,  recon¬ 
struction  time  can  possibly  be  reduced  by  grouping  coil  elements 
together  before  digitization  to  reduce  the  overall  number  of  data  chan¬ 
nels  [28,29].  However,  as  commercial  systems  continue  to  improve  in 
the  future,  reconstruction  times  are  expected  to  decrease  and  make 
this  less  of  an  issue. 

Another  possible  disadvantage  in  the  use  of  incremental  coils  locat¬ 
ed  somewhat  distantly  from  the  FOV  of  interest  is  that  motion  of  distant 
objects  might  alias  into  the  reconstructed  FOV  owing  to  the  high  sensi¬ 
tivity  of  the  incremental  coils  to  the  moving  object.  The  evaluation  per¬ 
formed  in  this  work  indicated  that  this  was  not  a  significant  problem,  if 
artifact  were  present  and  could  be  associated  with  some  specific  distant 
body  region,  then  conceivably  the  reconstruction  could  be  repeated 
with  data  from  the  coils  nearby  to  that  region  excluded. 

In  this  work  data  from  each  32-coil  DCE-MRI  acquisition  were  recon¬ 
structed  two  ways,  the  first  using  all  32  coils  and  the  second  using  only 
12  coils,  and  the  results  compared.  This  study  design  eliminated  the  ad¬ 
ditional  uncertainty  and  expense  associated  with  a  study  in  which  each 
subject  would  have  been  imaged  twice  with  contrast  material,  once 
with  each  coil  set.  We  validated  this  approach  by  comparing  12-coil  re¬ 
construction  from  32-  and  12-coil  acquisitions  and  showed  negligible 
difference  in  SNR  and  g-factor. 

This  study  evaluated  two  vendor-provided  configurations  of  the  coil 
arrays  available  with  the  MRI  system.  It  is  possible  that  other  configura¬ 
tions  might  provide  improved  performance.  Specifically,  acceleration 
applied  along  a  particular  direction  benefits  from  coil  elements  which 
tend  to  face  each  other  along  that  direction.  Electrically  combining  the 
two  most  lateral  elements  of  the  GEM  array  into  one  virtual  coil,  as 
done  by  the  vendor,  might  not  be  optimal  for  this  array  for  L/R  acceler¬ 
ation.  Also,  other  styles  of  receiver  coil  arrays  may  provide  improved 
performance,  such  as  arrays  placed  in  better  proximity  or  wrapped 
around  the  pelvis.  Incorporation  of  an  endorectal  coil  can  be  expected 
to  provide  some  improvement  in  SNR  near  that  coil.  In  this  work  al¬ 
though  an  endorectal  coil  was  placed  in  some  subjects  as  seen  in 
Fig.  7,  it  was  active  only  for  sequences  other  than  DCE-MRI. 

The  quantitative  analysis  distinguished  between  SNR  improvement 
solely  due  to  the  increase  in  the  number  of  coils  in  an  unaccelerated 
scan  (red  dots  of  Fig.  4)  as  well  as  that  due  to  improvement  in  g-factor 


Fig.  8.  (A)  Photograph  of  phantom  used  for  assessment  of  prostate  MRI.  Superior/inferior 
levels  of  bladder  and  prostate  inclusions  are  noted.  (B)  Axial  image  of  the  phantom  at  the 
level  of  prostate  inclusion  (dashed  circle)  showing  relative  percent  improvement  in  SNR  of 
32-  vs.  12-coil  acquisition.  (C)  Coronal  T2SE  images  of  the  phantom  acquired  using  12- 
(left)  and  32-coil  (right)  acquisitions.  The  observed  9.7%  SNR  increase  is  difficult  to 
visualize,  but  the  images  show  the  size  and  location  of  the  prostate  inclusion.  ROl  used 
for  SNR  measurement  is  indicated  as  dashed  ellipse  (left  image). 
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statistics  (increment  between  red  dot  and  black  triangle  for  each  pa¬ 
tient)  at  the  specific  acceleration  factor  used.  The  first  of  these  should 
benefit  any  sequence.  The  32-coil  array  provided  improvement  in 
both  SNR  and  in  g-factor  statistics  vs.  the  12-coil  array.  Both  are  impor¬ 
tant  in  acquisitions  such  as  DCE-MRI  in  which  acceleration  is  used.  It  is 
interesting  that  the  level  of  improvement  of  SNR  (red  circles  of  Fig.  4) 
appeared  to  correlate  positively  with  BMI. 

Extending  this  to  sequences  beyond  DCE-MRI,  the  phantom  study 
results  in  Table  4  indicate  that  both  T2SE  and  DWI  are  expected  to  ben¬ 
efit  with  improved  SNR  with  32-coil  acquisition.  Consistent  with  the  re¬ 
sults  shown,  sequences  which  are  particularly  signal-starved,  such  as 
high-b  value  DWI,  might  in  practice  benefit  more  from  the  SNR  increase. 
This  can  potentially  be  used  to  reduce  the  level  of  averaging  commonly 
done  in  DWI. 

The  benefit  of  increased  SNR  can  be  used  for  DCE-MRI  in  various 
ways.  As  stated,  we  initially  investigated  but  opted  not  to  go  to  faster 
frame  rates.  We  have  developed  a  version  for  assessment  of  biochemical 
recurrence  with  higher  spatial  resolution  (0.76  x  0.87  x  2.24  mm3  vs. 
the  0.86  x  1.15  x  3.0  mm3  of  this  work)  with  a  15  s  vs.  6.6  s  frame  time. 

In  addition  to  the  above-described  extended  reconstruction  time 
and  possible  motion  artifact  from  distant  objects,  another  limitation  of 
this  work  is  that  further  study  is  necessary  to  determine  whether  pros¬ 
tate  lesion  detection  and  characterization  are  improved  with  32-coil 
DCE-MRI.  Also,  study  of  the  impact  in  T2SE  and  DWI  sequences  should 
be  extended  to  a  detailed  patient  series  for  a  more  complete  assessment. 

In  summary,  for  the  same  level  of  acceleration,  the  acquisition  and 
reconstruction  of  data  from  32  receiver  coils  encompassing  the  pelvis 
provides  superior  measured  and  perceived  SNR  in  dynamic-contrast- 
enhanced  prostate  MRI  vs.  use  of  12  receiver  coils. 

Supplementary  data  to  this  article  can  be  found  online  at  http://dx. 
doi.org/1 0.1 01 6/j.mri.201 7.01 .01 7. 
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